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We study the magnetically-induced phonon splitting in cubic ACr204 (A=Mg, Zn, Cd, Hg) 
spinels from first principles, and demonstrate that the sign of the splitting, which is experimentally 
observed to be opposite in CdCr204 compared to ZnCr2 04 and MgCr2 04, is determined solely by 
the particular magnetic ordering pattern observed in these compounds. We further show that this 
interaction between magnetism and phonon frequencies can be fully described by the previously 
proposed spin-phonon coupling model that includes only the nearest neighbor exchange. Finally, 
using this model with materials specific parameters calculated from first principles, we provide 
additional insights into the physics of spin-phonon coupling in this intriguing family of compounds. 


I. INTRODUCTION 

The interplay of spin and lattice degrees of freedom 
can lead to a variety of fundamentally and technologically 
interesting phenomena including the spin-Jahn-Teller ef¬ 
fect in frustrated magnets , 1,2 magnetocapacitance , 3 and 
the linear magnetoelectric effect . 4,5 One signature of this 
interplay is the influence of magnetic order on the vibra¬ 
tional spectrum of a material. In many transition metal 
oxides the spin correlations shift the phonon frequen¬ 
cies, and lead to the so called magnetodielectric effect . 6,7 
Furthermore, if the long-range antiferromagnetic (AFM) 
order reduces the crystal symmetry, the onset of anti¬ 
ferromagnetism can result in a substantial splitting of 
phonon frequencies that are degenerate in the paramag¬ 
netic (PM) phase, even when the change in the crystal 
structure is undetectable . 8-12 This phonon anisotropy is 
a non-relativistic effect which originates from the changes 
in hybridization due to spin ordering. In particular, the 
phonon splitting can be phenomenologically explained by 
a dependence of the exchange interactions on the atomic 
positions . 13,14 

Chromium spinels ACr204 (A=Mg, Zn, Cd) are a par¬ 
ticularly interesting class of frustrated antiferromagnets 
that exhibit strong spin-phonon coupling. In the PM 
phase, group theory predicts, and experiments confirm, 
the presence of four triply degenerate infrared (IR)-active 
phonon modes. Below the Neel temperature, however, 
one of these phonon modes undergoes a large splitting 
into a singlet and a doublet 10,15-17 . This feature and 
its magnitude in was argued to be a consequence of a 
dominant role of the nearest-neighbor (nn) direct Cr- 
Cr exchange interaction 10 . Fennie and Rabe developed 
in Ref. [ 11 ] a general approach to incorporate mate¬ 
rial specific information from first principles into spin- 
phonon coupling models. They demonstrated that the 
spin-phonon coupling model with only nearest-neighbor 
(nn) exchange interaction and parameters derived from 
first principles provide a full description of experimen¬ 
tally observed magnetically-induced phonon splitting in 


ZnC^C ^. 11 This model has later been successfully ap¬ 
plied to many other systems as well . 18 21 

Magnetically-induced phonon splitting have been also 
observed in other Cr-based spinels including MgCr 2 04 
and CdCr204. Interestingly, the sign of the phonon 
slipping observed for MgCr204 and ZnC^C^ (kinglet > 
^doublet) 10,17 is opposite to that observed in CdCr204 
(^doublet > ^singlet)- 15,16 In all three of these compounds, 
the sign of the nn exchange interaction is the same, how¬ 
ever its magnitude compared to further neighbour ex¬ 
changes is dramatically different. In the Mg and Zn 
compounds the nn exchange interactions are two orders 
of magnitude larger than all other exchange interactions, 
while in CdC^C^ the nn interaction is the same order 
as the second neighbor interaction. Based on this fact 
Kant et al 17 concluded that the spin-phonon coupling 
model with only nn exchange interaction cannot explain 
the magnetically induced phonon anisotropy in AC^CU 
spinels, and instead proposed that the phonon splitting 
is generally controlled by a nondominant, next neighbour 
exchange interaction. 

In this paper we use first principles calculations to 
study the magnetically induced phonon anisotropy of the 
zone-center polar modes in AC^CU (A=Mg, Zn, Cd, 
Hg) spinels. We show that the different magnetic or¬ 
derings characteristic for these spinels lead to different 
signs of the phonon splitting. In particular, we explain 
the opposite sign observed for ZnCr204 and MgCr204 
compared with CdCr 2 04 which have distinct magnetic 
ground states. We find that the spin-phonon coupling 
model of Ref. [ 11 ] with only the nn exchange interactions 
can very accurately describe ab initio values of phonon 
frequencies for all the spinel compounds we considered. 


II. METHODS 

The first principles calculations were performed using 
the density functional theory (DFT) within the rotation- 
ally invariant DFT+U method 22 and the PBEsol approx¬ 
imation to the exchange-correlation functional 23 . Simi- 
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FIG. 1. (a) Cubic crystal structure of ACr 2 04 spinels con¬ 
sisting of Cr-centered octahedra and A-centered tetrahedra. 
(b) Magnetic exchange couplings up to the third Cr neigh¬ 
bors; note that there are two nonequivalent types of third 
neighbors which have distinct exchange parameters: J 3 and 
J 3 . Brown, blue and red spheres denote Cr, A and O atoms, 
respectively. 


larly as in Ref. 11 we used U = 3 eV and J = 0.9 eV, 
the parameters that accurately reproduce photoemission 
spectra and band gaps in sulfur Cr spinels 24 . 

The Kohn-Sham equations were solved using the pro¬ 
jector augmented wave method 25 as implemented in the 
VASP code 26,27 with the plane wave cutoff of 500 eV 
and 6x6x6 T-centered k-point mesh in the primitive 
unit cell for the cubic structure (for larger cells the k- 
point mesh was scaled accordingly). Phonon frequen¬ 
cies and eigendisplacements were calculated using the 
frozen phonons method using symmetry adapted modes 
obtained from the ISOTROPY package 28 . The spin-orbit 
coupling was neglected in the calculations. Structural re¬ 
laxations were performed in the ferromagnetic (FM) state 
that preserves the cubic symmetry. 


III. CRYSTAL AND MAGNETIC STRUCTURE 

At high temperatures A&2O4 spinels have a cubic 
(Fd 3 m) structure where A 2+ ions are in tetrahedral oxy¬ 
gen environment and form the diamond lattice, while 

Cr 3 + 

ions are surrounded by octahedral oxygen cages 
and form a pyrochlore lattice, see Fig. la. The calcu¬ 
lated structural parameters are shown in supplementary 
materials 29 . 

The octahedral crystal field splits the Cr 3 d orbitals 
into a lower-lying t^g triplet and a higher-energy e g dou¬ 
blet. Cr 3+ has three outer electrons that fill the major¬ 
ity t2 g states which results in a net Cr spin 5 = 3 / 2 . We 
found the exchange interaction parameters between Cr 
spins by fitting ab initio energies of different collinear 
magnetic configurations to the Heisenberg Hamiltonian 


A 

Ji Ji J 3 J3 


JahM m 

12 3 4 

12 3 4 

Mg 

3.81 -0.07 0.05 0.18 

4.5 2.4 0.4 0.3 

-0.1 -0.1 -0.7 -0.7 

Zn 

3.81 -0.08 0.14 0.16 

2.0 5.0 0.8 0.2 

-0.0 -0.1 -0.8 -0.7 

Cd 

0.29 -0.10 0.08 0.15 

1.3 3.3 0.7 0.7 

-0.0 -0.0 -0.5 -0.2 

Hg 

-0.58 -0.00 0.21 0.13 

0.8 3.2 0.9 0.6 

-0.0 -0.0 -0.5 -0.1 


TABLE I. The exchange couplings (in meV) and L7 a _l, || / M 
parameters (in cm -1 ) calculated for different ACr 204 spinels. 


h = J2 JijSi • s, (1) 

where the summation is over Cr ions, Jij are the ex¬ 
change parameters, and S i is the unit vector indicating 
the direction of the spin at Cr site i. We considered ex¬ 
change parameters up to the third neighbors (see Fig. 
lb) as further neighbors are known to have negligible ex¬ 
change couplings 30 . The results are presented in Table 
I. 

The nn exchange parameter, Ji, is a dominant inter¬ 
action for all compounds. This coupling arises from the 
competition between AFM direct exchange and FM 90 ° 
superexchange 31 . For A ions with small ionic radii, like 
Mg or Zn, the direct exchange mechanism dominates 
resulting in a strong AFM J\. However, for larger A 
ions the lattice parameter and the nn Cr-Cr distance 
increases 29 which diminishes the direct exchange con¬ 
tribution. In particular, for CdCr 2 04 the AFM J\ is 
reduced by an order of magnitude while for HgCr 2 04 
the superexchange contribution overcomes the direct ex¬ 
change resulting in a (small) FM J\. 

Exchange couplings beyond nn originate from higher- 
order superexchange processes 32 and, in general, are 
smaller than J\. However, while for MgCr 2 04 and 
ZnCr 2 04 these interactions are negligible compared to 
the nn exchange, for CdCr 2 04 and HgCr 2 C>4 compounds 
the J3 and J3 exchanges become relevant. 

The AFM nn exchange interaction is frustrated on the 
pyrochlore lattice since spins forming a tetrahedron can¬ 
not be all antiparallel to each other. The energy due to J\ 
is minimized when for all tetrahedra the total spin is zero, 
i.e. in each tetrahedron two spins are parallel while the 
other two point in the opposite direction. There are, how¬ 
ever, many such two-up-two-down configurations which 
can be different in different tetrahedra leading to infinite 
degeneracy. Consequently, the magnetic ground state 
is determined by further exchange couplings 33 , magne- 
toelastic effects 1,2 or relativistic interactions 34 leading to 
complicated, often noncollinear, orderings. 

There are two primary collinear magnetic orders that 
are relevant for ACr 2 04 spinels. These are shown in Fig. 
2 and we denote them as AFM-I and AFM-II. The AFM- 
I ordering is similar to the true magnetic ground state in 
ZnCr 2 04 2 and MgCr 2 04 35 while the AFM-II approxi- 
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FIG. 2. Two collinear magnetic orderings relevant for ACr 2 04 
spinels, (a) AFM-I order that is similar to the true magnetic 
ground state in ZnCr 204 and MgCr 2 04 - (b) AFM-II order 
that approximates the spin order in CdCr 2 04 - Brown and 
red spheres denote Cr and O atoms, respectively. 

mates the spin order in CdCr 2 C>4 33 . Both spin orderings 
satisfy the two-up-two-down rule in each tetrahedron but 
they differ in relative orientations of spins in neighboring 
tetrahedra. In particular, for AFM-I the nearest neigh¬ 
bors in the xy plane are parallel while for AFM-II they 
are antiparallel. As we will see below, this difference has 
a profound effect on the magnetically induced phonon 
anisotropy. 


IV. PHONON FREQUENCIES 

We now focus on the influence of magnetic order on 
the zone-center polar phonons. In addition to the AFM- 
I and AFM-II orderings which are relevant for this class 
of compounds (see above), we also considered the FM 
order since which has the same cubic symmetry as the 
PM state. 

In the FM state there are four triple degenerate polar 
phonon modes, each transforming according to the T\ u 
irreducible representation of the Oh cubic point group. In 
order to compute the corresponding phonon frequencies 
we considered symmetry-adapted T\ u modes, / n?a . Here 
n = 1,2, 3 , 4 , 5 is the mode number (in addition to four 
polar modes we need to include the acoustic mode that 
also has a T\ u symmetry) and a = x,y,z labels the row 
of Tiu such that f n , a transforms as a vector along the 
a axis. The symmetry-adapted modes, /i ?a , / 2?cn and 
/4 ?cn involve atomic displacements along a of the entire 
A, Cr, and O sublattice, respectively. On the other hand, 
atomic displacements associated with and /5 ?a take 
place in the plane perpendicular to the a axis and in¬ 
volve chromium and oxygen atoms, respectively (see the 
supplementary materials for more details). 

Condensing the symmetry-adapted modes for a given 
a and evaluating the Heilman-Feynman forces for the FM 
state, we constructed the 5 x 5 dynamical matrix block. 
Matrix diagonalization leads then to four nonzero phonon 
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FIG. 3. Calculated T lu phonon frequencies (in cm -1 ) of dif¬ 
ferent ACr 2 04 spinels for different magnetic orderings. In the 
FM state (middle) there are four optical T\ u modes and each 
one is three-fold degenerate. In the AFM-I (left) and AFM-II 
(right) states each triplet splits into a singlet and a doublet. 
For the AFM-I (AFM-II) ordering the singlet (doublet) has 
higher frequency for all compounds and for all phonon modes. 


frequencies corresponding to the four polar phonon 
modes, see Fig. 3 (middle). As expected, the frequencies 
are independent of a leading to three-fold degeneracy of 
each mode. 

For the AFM-I and AFM-II orderings the symmetry 
is lowered to tetragonal D^h and D4 point groups, re¬ 
spectively, with the tetragonal direction chosen to be 
along the z axis. In both cases the T\ u representation 
becomes reducible resulting in a splitting of the triple 
degenerate polar phonon modes according to T\ u -A 
A 2u ® E u for AFM-I and T\ u for AFM-II. The 

one-dimensional A 2u and A 2 irreducible representations 
transform as a vector along the z axis, while the two- 
dimensional E u and E irreducible representations trans¬ 
form as a vector in the xy plane. 

In order to evaluate this splitting we calculated the po¬ 
lar phonons for the AFM-I and AFM-II orderings by di¬ 
agonalizing the three 5 x 5 blocks of the dynamical matrix 
using the Tiu symmetry-adapted modes, / n?Q , 36 . Since 
we are interested in the splitting generated by the spin 
pattern alone, similarly as in Ref. 8 we neglected the 
magnetically induced tetragonal distortion of the crystal 
and performed calculations for the cubic structure found 
using FM configuration. 

The results are shown in Fig. 3 . The phonon fre¬ 
quencies obtained from the a = x and a = y dynamical 
matrix blocks are equal and form E u or E doublets. On 
the other hand, the phonon frequencies obtained from 
the a = z block are different and correspond to A 2u 
and A 2 singlets. The phonon splitting is the largest for 
the second lowest frequency mode (except for MgCr 2 04 
where it is the lowest frequency mode that has the high- 
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est splitting). In particular, for ZnCr 2 04 in AFM-I state 
it becomes as large as 41 cm -1 . Interestingly, the mag¬ 
nitude of the splitting for the AFM-I order is, in general, 
about twice larger than for the AFM-II state. 

The most important feature, however, is that for all 
considered compounds and for all the modes the singlet 
has a higher energy for the AFM-I state while for the 
AFM-II configuration it is the doublet that has a higher 
energy. This is in agreement with the experimentally ob¬ 
served sign reversal of the phonon splitting for ZnCr 2 04 
and MgCr 2 04 compared with CdCr 2 04 10,15 17 since the 
former ones have a ground state similar to AFM-I 2 while 
the magnetic ordering of the latter can be approximated 
by AFM-II 33 . 

V. SPIN-PHONON COUPLING MODEL 

In order to better understand the results of our first 
principles calculations we employ the spin-phonon cou¬ 
pling model in which the Tin block of the force-constant 
matrix for an arbitrary magnetic state is given by 11 

__ 02 t 

c™,** = c-™ + 4 £ y (Si • s,-) (2) 

j O'Jna^'Jn'a' 

Here, (S i • S j) is the spin correlation function that is 
1 for the FM ordering and either 1 or -1 for the AFM 
ordering. < 7 ™, is the force-constants matrix in the PM 
phase. Note that the latter has a Oh cubic symmetry 
and thus it doesn’t depend on the T\ u row indices a and 
a'. 

The above expression can be further simplified by us¬ 
ing the symmetry of the magnetic state. In particu¬ 
lar, the AFM-I and AFM-II (as well as FM) orderings 
don’t induce couplings between different rows of T\ u so 
the force-constant matrix is diagonal in the row indices: 
Cna.n'oi' = C na ^ a 5 a ^ = c n>n /(a). Consequently, we 
only need d 2 Jij/df na df n ' a . Considering only the nn 
exchange interaction, there are only two types of such 
derivatives: 11 J" n ,_ L = d 2 Jij/df na df n ' a V f i<7 - • a = 0 and 
J nn'\\ ~ < 9 2 Jij/dfnadfn'a V Tij * & A 0 , where Yij is the 
vector linking nn sites i and j and a is the vector along 
the a axis. Therefore, we can write 


(7 n ,„4a)=C™, +4J" n , i ^(S i 

r_L 

• s,) 


+ 4 ^nn'|| N( S » 

•Si) 

(3) 


where the first summation is over the two nn in the 
plane perpendicular to <a, and the second summation is 
over the other four nn. For the three magnetic orderings 
considered, we obtain 

C™(a = y, z) = C™ + 8J" n , ± + 16 J" n ,|| ( 4 ) 


= x,y)= C™ - 8 J"„, ± ( 5 ) 


= z) = C™ + 8J'n n i± - 16 J" nl (6) 


C™(a = x,y) = C™ - S./",, r ( 7 ) 


Ct™- n (a = ^) = C™-8J" n , ± (8) 

The above equations explicitly demonstrate that in the 
FM state we have a three-fold degeneracy with respect to 
a and that in the AFM-I and AFM-II states these triplets 
split into a doublet (a = x, y ) and a singlet (a = z). 

The parameters ( 7 ™, and J" n , y were fitted to 

the ab initio force-constant matrices evaluated for FM, 
AFM-I, and AFM-II orderings. Essentially perfect fitting 
was obtained with the misfit lower than 0.03 meV/A and 
corresponding phonon frequencies within 1 cm -1 of first 
principles values (see Supplementary Materials 29 ). This 
indicates that the spin-phonon coupling model with only 
nn exchange coupling provides an excellent description 
of the effect of magnetic order on phonon frequencies in 
ACr 2 04 spinels. 

Explicit forms of C PM , J", and J" force-constant ma¬ 
trices for different ACr 2 04 spinels are shown in Supple¬ 
mentary Materials 29 . We find that for all compounds 
e/ggj_ is positive and significantly larger than any other 
element of and Jy' matrices. As discussed in Ref. 
11, the anomalously large value of J^' 3± originates from 
the exponential form of the direct exchange contribu¬ 
tion (Jd) to J\. Indeed, Jd = Ae~ aDcr ~ Cr where A 
and a are positive constants and Dc r -Cr is the nn Cr- 
Cr bond length. The only partner function that signif¬ 
icantly affects Dcr-Cr is /3a (with a perpendicular to 
the bond) 10 resulting in large J 33 _^. This explanation is 
consistent with the fact that the J 33± element is similar 
for MgCr 2 04 and ZnCr 2 04 compounds but it decreases 
with the size of A ion due to diminished role of the direct 
exchange mechanism. Note also that the positive sign of 
J 33± is the direct consequence of the exponential depen¬ 
dence of Jd on the atomic displacements which requires 
the second derivative to have the same sign as J^. 

The second largest element among and matri¬ 
ces is J 3 hi_ ( or J$3±)- Since f^ a corresponds to the dis¬ 
placement of O sublattice which modifies Cr-O-Cr an¬ 
gle, this shows that the superexchange mechanism also 
contributes to the spin-phonon coupling. However, the 
superexchange contribution to the spin-phonon coupling 
is always significantly smaller than the direct exchange 
contribution. This remains true even for HgCr 2 04 com¬ 
pound where the superexchange interaction is stronger 
than the direct exchange coupling. We believe that this 
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relative ineffectiveness of the superexchange in generat¬ 
ing a spin-phonon coupling is a generic feature and is 
due to the fact that this mechanism doesn’t depend so 
strongly on atomic displacement as the direct exchange 
mechanism. This feature plays also an important role in 
the success of our spin-phonon coupling model where we 
considered only nn exchange interaction. Indeed, while 
for for MgCr204 and ZnC^CU J\ is at least order of 
magnitude larger than other exchange parameters, in 
the case of CdC^CU and HgC^CU compounds the J 3 
and J3 couplings are not negligible and our approxima¬ 
tion works only because these couplings originates from 
superexchange processes and have weak dependence on 
atomic displacements. 

Having established and understood the validity of our 
spin-phonon coupling we can now use it to provide an ad¬ 
ditional insight into magnetically induced phonon split¬ 
ting. Since J" n , ± and J” n , y are much smaller that the 
elements of the paramagnetic force-const ant matrix, the 
phonon frequencies can be written as 


ti\(a) k w™ + 

|m I ^X E< S * • S ;> + ^All E< s * • S ;> ) 

A \ r,| / 


( 9 ) 


Here, cj™ is the paramagnetic phonon frequency and 
for each phonon mode we introduced: y = u\j " ^u\ 

where u\ are the paramagnetic dynamical matrix eigen¬ 
vectors and J" || is the dynamical matrix correspond¬ 
ing to J” n rj_ || force-constant matrix. The magnetically- 
induced phonon splittings for the two principle AFM or¬ 
ders are then given by 


A^™- 1 EE CD™ 1 ! 


a = z) — cjy 


(a = x,y) 




AFM-II _ 


a = z - cm 


l (a = x,y) 


-4JTxK M +4^,M M 


The ratio y/u;™ characterize the strength of the 
magnetic contribution to the phonon frequencies. Table 
I shows these parameters for different AC^CE spinels. 
We can immediately observe that the only appreciable 
ratio are for the two lowest-frequency modes 

(A =1,2) and both are always positive. It follows from 
Eqs. ( 10 ) and ( 11 ) that the phonon splitting is the largest 
for the two lowest-frequency modes and it is positive 
(negative) for the AFM-I (AFM-II) orderings. This fea¬ 
ture is a direct consequence of J 33± being positive and 
dominant among other elements of J" and J" since the 
A =1,2 modes have the largest content of the f 3a part¬ 
ner function (see Supplementary Materials 29 ). Therefore, 


according to the discussion above, different signs of the 
phonon splittings that we found from first principles for 
the lowest phonon modes are ultimately related to the 
dominant role of the direct exchange mechanism in gen¬ 
erating the spin-phonon coupling. 

Interestingly, for A = 3,4 are also always pos¬ 

itive while J’x\\/ UJ \ M for all modes are consistently neg¬ 
ative. This results in the sign of the phonon splitting 
to be positive (negative) for the AFM-I (AFM-II) states 
for the two-highest phonon modes as well. These fea¬ 
tures, however, is difficult to explain microscopically due 
to small values of the splittings. In fact, the spliting of 
the two-highest phonon modes is too small to be seen in 
experiments. 

According to Eqs. ( 10 ) and (11) we can write 

Aw™ 1 « - Aw™ (12) 

Therefore, the phonon splittings for AFM-I and AFM- 
II orders are always opposite and the latter is approxi¬ 
mately half of the former. This is a general result which 
follows directly from the validity of the nn spin-phonon 
coupling model and it is independent on the signs and 
sizes of the exchange second derivatives. Note that the 
above relation is well satisfied by first principles data 
(Fig. 3 ) which again demonstrates validity of the model. 
Note, however, that this relation doesn’t tell us for which 
ordering the phonon splitting is positive. In order to an¬ 
swer this question more microscopic analysis (as above) 
is needed. 


VI. CONCLUSIONS 

In summary, we investigated the effect of magnetic 
ordering on phonon frequencies of ACr204 spinels us¬ 
ing first prinicples electronic structure calculations. We 
found that our ab initio results are very well described by 
the spin-phonon coupling model with only nn exchange 
coupling 11 . Both the model and first principles calcu¬ 
lations show that a specific type of spin ordering has 
a crucial effect on magnetically induced phonon split¬ 
ting. In particular, we found that the different magnetic 
states observed in different spinels lead to the opposite 
signs of the phonon splittings observed in ZnC^CU and 
MgCr2 04 compounds compared to CdC^C^. This fea¬ 
ture is a result of an important role played by the direct 
exchange mechanism in generating the spin-phonon cou¬ 
pling in these materials. 
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I. STRUCTURAL DATA 


Compound 

a 0 (A) x 

Dcr-Cr (A) 

a (°) 

MgCr20 4 

8.372 0.262 

2.960 

95.4 

ZnCr20 4 

8.351 0.262 

2.952 

95.4 

CdCr 2 0 4 

8.628 0.269 

3.051 

98.8 

HgCr20 4 

8.703 0.271 

3.077 

99.8 


TABLE I. Structural parameters of different AC^CL spinels 
in the Fd3m cubic phase calculated for the ferromagnetic spin 
configuration. In addition to the lattice parameter (ao) and 
the fractional coordinate of oxygen atoms (x) we also show the 
nearest-neighbor Cr-Cr distance ( Dcr-Cr ) and the Cr-O-Cr 
angle (a). 
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ZnCro0 4 


C PM - 


9.269 

-2.470 

-4.047 

-2.888 

-3.815 

-2.471 

20.804 

4.428 

-13.475 

-1.104 

-4.049 

4.426 

31.416 

-1.106 

-10.057 

-2.888 

-13.467 

-1.108 

10.967 

2.690 

-3.817 

-1.105 

-10.049 

2.690 

16.631 


II. SYMMETRY ADAPTED MODES 


The symmetry-adapted modes f UjZ (n = 1,2,3,4,5) 
are illustrated in Fig. SI. The f n , x and f n ^ y modes can 
be obtained from f n , z by the 90° rotation around the y 
and the x axis, respectively. 
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III. FITTING RESULTS 

Fitted values of the parameters of the spin-phonon cou¬ 
pling model for different ACr 2 0 4 compounds are shown 
below. The units are meV/A. 

MgCr 2 0 4 
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FIG. SI. Symmetry-adapted modes f n , z (n = 1,2, 3,4, 5): (a) /i ;Z , (b) f 2 , z , (c) f 3 ,z, (d) f^z, (e) /s,*. Blue, brown, and red 
spheres denote A, Cr, and O atoms, respectively while the arrows indicate atomic displacements. Symmetry-adapted modes 
fn,x ( fn, y ) can be obtained by the 90° rotation around the y (x) axis. 
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The phonon frequencies for FM, AFM-I, and AFM-II 
orderings evaluated from the spin-phonon coupling model 
are shown in Table II. As expected, the results agree very 
well with first principles calculations (see Fig. 3 of the 
main manuscript). In addition, we show the paramag¬ 
netic phonon frequencies. 
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Compound 

A 


FM 

AFM-I 

AFM-II PM 

x,y,z 

x,y 

z 

x,y 


x,y,z 

MgCr2 04 

1 

261 

224 

261 

245 

224 

245 

2 

424 

406 

425 

415 

406 

415 

3 

477 

479 

488 

484 

479 

481 

4 

612 

615 

623 

619 

615 

616 

ZnCr204 

1 

189 

173 

189 

182 

173 

182 

2 

379 

340 

381 

360 

340 

360 

3 

499 

499 

511 

505 

499 

502 

4 

599 

603 

611 

607 

603 

604 

CdCr 2 0 4 

1 

151 

141 

151 

146 

141 

146 

2 

369 

343 

370 

356 

343 

356 

3 

473 

471 

481 

476 

471 

474 

4 

598 

594 

601 

597 

594 

596 

HgCr204 

1 

108 

101 

108 

105 

101 

105 

2 

357 

332 

358 

345 

332 

345 

3 

475 

471 

483 

477 

471 

475 

4 

575 

571 

577 

574 

571 

573 


TABLE II. Phonon frequencies for different magnetic order¬ 
ings calculated from the spin-phonon coupling model. The 
units are cm -1 . 


IV. PARAMAGNETIC DYNAMICAL MATRIX 
EIGENVECTORS 


Compound 

A 

n 

1 

2 

3 

4 

5 

MgCr 2 04 

1 

0.38 

0.05 

0.36 

0.01 

0.20 

2 

0.45 

0.06 

0.45 

0.01 

0.03 

3 

0.00 

0.27 

0.07 

0.42 

0.24 

4 

0.04 

0.08 

0.12 

0.23 

0.52 

ZnCr204 

1 

0.58 

0.16 

0.12 

0.07 

0.07 

2 

0.13 

0.04 

0.65 

0.01 

0.16 

3 

0.00 

0.28 

0.06 

0.47 

0.19 

4 

0.01 

0.07 

0.16 

0.18 

0.58 

CdCr 2 0 4 

1 

0.51 

0.24 

0.08 

0.11 

0.06 

2 

0.07 

0.03 

0.78 

0.02 

0.10 

7 

0.00 

0.30 

0.03 

0.44 

0.22 

4 

0.01 

0.06 

0.11 

0.20 

0.62 

HgCr204 

1 

0.42 

0.32 

0.05 

0.18 

0.03 

2 

0.03 

0.02 

0.81 

0.02 

0.12 

7 

0.00 

0.33 

0.02 

0.48 

0.16 

4 

0.00 

0.05 

0.12 

0.14 

0.69 


TABLE III. Squared projections of the paramagnetic dynami¬ 
cal matrix eigenvectors, u\, on the symmetry adapted partner 
functions, f na . 



